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A study by W. R. Magro and D. M. Ceperley [Phys. Rev. Lett. 73, 826 (1994)] has shown 
that the ground state of the two-dimensional fluid of charged bosons with logarithmic interactions 
is not Bose-condensed, but exhibits algebraic off-diagonal order in the single-particle density matrix 
p(r). We use a hydrodynamic Hamiltonian expressed in terms of density and phase operators, in 
combination with an /-sum rule on the superfluid fraction, to reproduce these results and to extend 
the evaluation of the density matrix to finite temperature T. This approach aUows us to treat the 
hquid as a superfluid in the absence of a condensate. We find that (i) the off-diagonal order arises 
from the correlations between phase fluctuations; and (ii) the exponent in the power-law decay of 
p(r) is determined by the superfluid density UsiT). We also find that the plasmon gap in the single- 
particle energy spectrum at long wavelengths decreases with increasing T and closes at the critical 
temperature for the onset of superfluidity. 
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I. INTRODUCTION 

The fluid of point charges interacting via a logarithmic Coulomb potential in a strictly two-dimensional (2D) 
space is a basic model in statistical mechanics with main relevance to the theory of vortex fluctuations in superfluid or 
superconducting filmsS. Magro and Ceperleyll have discussed the bosonic ground state of this model system. They first 
used a sum-rule argument from the work of Pitaevskii and Stringarii to show that the zero-point fluctuations associated 
with long-wavelength plasmons rule out the presence of a Bose-Einstein condensate even at zero temperature. They 
proceeded to evaluate the single-particle density matrix p{r) by quantum Monte Carlo (QMC) methods and showed 
that the ground state exhibits algebraic off-diagonal long-range order: the decay of p{r) with increasing distance r is 
through the power law r~'^'/'^, where is a coupling strength parameter determined by the areal density n. 

As Magro and Ceperleyi emphasize, although the ground state of the 2D Bose Coulomb fluid (2D-BCF) is not 
condensed, superfluidity is nonetheless possible. If the model is indeed superfluid, it may then provide an ideal system 
in which to study the differences between superfluidity and Bose-Einstein condensation. 

A number of structural and collective dynamical properties have been calculated for the 2D-BCF within approximate 
theoretical schemed, but no other studies seem to have been made of its momentum distribution and one-body density 
matrix beyond that of Magro and Ceperley. In the present paper we evaluate these properties for the charged fluid in 
the weak-coupling regime corresponding to high density, by adopting a hydrodynamic reduction of its Hamiltonian as 
previously proposed by Popov! for 2D neutral Bose fluids (see also Men^) . A crucial point of our approach is to allow 
for a difference between the long-wavelength dispersion relations of single-particle and collective excitations at finite 
temperature T, through the use of distinct /-sum rules on the superfluid and the total particle number densitiesB. 
We recover by this approach the results of Magro and Ceperleyi for the power-law decay of p{r) in the ground state 
and obtain their extension to finite temperature. We show that within this approach the 2D-BCF is indeed superfiuid 
at T = and that in the weak-coupling regime a slowly declining superfiuid fraction persists as the temperature is 
raised. In fact, the superfluid fraction enters to determine the value of the exponent for the power-law decay of p{r). 
The presence of the superfluid density in the power-law decay of the density matrix in a neutral 2D Bose gas was 
already argued by Popovi. 

Our present progress in understanding the properties of the 2D-BCF borrows from the theory of a quasicondensate 
state in the neutral 2D Bose gas, that was developed by a number of authorsi@~0 and has also been used within 
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a Bogoliubov approach to describe the 2D fluid of charged bosons with e^/r interactionala. In essence, the local 
properties of a quasicondensate, over distances that are shorter than the phase correlation length, are the same 
as those of a genuine condensate. However, in our treatment of the 2D-BCF we dispense with the notion of a 
quasicondensate fraction and base our arguments on the /-sum rule for the superfluid fraction. 

The paper is organized as follows. In Section II, after introducing the Hamiltonian of the charged fluid with its 
neutralizing background and recalling the /-sum rules for a superfluid from the work of Hohenberg and MartinB^ we 
present the hydrodynamic Hamiltonian as expressed in terms of density and phase operators. Following MengB we 
diagonalize it in the momentum representation and obtain a dispersion relation for single-particle excitations which 
has the Bogoliubov form but contains the superfluid density in place of the condensate density. In Section III we use 
the hydrodynamic Hamiltonian to determine the power-law decay of the one-body density matrix from the behavior 
of the correlations between phase fluctuations entering the single-particle Green's function. In Section IV we evaluate 
p{r) by using the resummation method adopted by other authorsi'i to deal with a fluid with a quasicondensate, and 
show that the power-law decay obtained from the single-particle Green's function is recovered. We then report in 
Section V numerical results for the superfluid density from Landau's formula and for the one-body density matrix 
as functions of temperature at various values of the coupling strength in the weak-coupling regime. Section VI 
concludes the paper with a brief summary. In an Appendix we show that the standard Bogoliubov approachEl would 
indeed yield that the condensate fraction is zero in the ground state of the 2D-BCF. 



II. HYDRODYNAMIC HAMILTONIAN AND SINGLE-PARTICLE EXCITATIONS 



The 2D-BCF is described by the Hamiltonian 

H = j drV^(r) - A* V^r) + ^ j dr j dY'^\v)^\v')V{\Y ~v'\)iI){y')iI){y) (1) 

(having set h = 1), where il^ir) is the field operator, /i is the chemical potential and V{r) is the interaction potential. 
This is the solution of the 2D Poisson equation V^F(r) — 27re^5(r), yielding V{r) = — e'^ln(r//o) where Iq is a 
reference length that we shall take as Iq = (me^)~^/^. The Fourier transform of the potential is Vk — 27re^/fc^, and 
in the following we shall set 14;=o = on account of a uniform background ensuring global charge neutralityQ. The 
coupling strength is measured by the dimcnsionless parameter r^, which is defined by = (2me^/7rn)^/^. 

A special role will be played in the following by a sum rule involving the superfluid density . We first recall that 
the usual /-sum rule0 involving the particle density n can be recast through the continuity equation into a sum rule 
on the longitudinal current-current response function xjj(k, w), 

duj Imxjj(k,^) ^ n_ 

For a charged fluid this is equivalent to the well-known plasmon sum rule on the longitudinal electrical current density. 
In addition to Eq. (2), a sum rule involving the response function of the superfluid velocity Vs holds in a superfluidElEj, 
which reads 

^.^^ d^Imx....(k,c.)^^_ (3) 

The role of this relation as a long-wavelength /-sum rule on the superfluid response becomes evident from its similarity 
to Eq. (2) when it is rewritten as a sum rule for the superfluid particle current density Js — ti^Vs. 

We are now ready to proceed to a hydrodynamic reduction of the Hamiltonian (1), following the work of PopovS 
and MenJ^ on neutral superfluids. This is obtained by making the transformation 



V'(r) = VpW exp [i$(r)], V^(i-) exp [-i$(r)] (4) 
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for the components of the field operators which correspond to wave number k below a cut-ofi^ kjo. Here p{r) and $(r) 
are the particle density and phase operators. The transformed Hamiltonian can be brought to a quadratic form by 
considering only small fluctuations in p{r) around a constant value po, 

p{r) = po + T]{r). (5) 

A term linear in 77(r) can be dropped by setting p = Vk=o = Oi and we obtain the hydrodynamic Hamiltonian 

Hh = J dr (8mn,)-i (V?7(r)) V (n,/2m) (V$(r))^ + ^ J dr J dr't/(|r - r'|)7/(r)77(r'). (6) 

Here, the constant po has been taken equal to ris in order to satisfy the sum rule (3), as we shall explicitly demonstrate 
at the end of this Section. As usual, the superfluid velocity is ~ V$(r)/m. We stress that the Hamiltonian in 
Eq. (6) describes only long-wavelength fluctuations of small amplitude in the superfluid density and velocity field. 
The Hamiltonian (6) can be diagonalizcd by first expanding the density fluctuation and phase operators in the form 

77(r) = ^ gkexp(ik-r), $(r) = ^ p_k exp (ik • r), (7) 

k(fc<fco) k(fe<feo) 



which yields 



Hh= ^ I (fc^/Smris) + <7kg-k + (nsA:^/2m)pkP-k 1 



k{k<ko) 



(8) 



We next set gk = (2Q;k) ^^^(ok + ^Lk) ^^^"^ P-k = — *ak(2ak) ^^^(^k — a^k) with [ak,ajj.,] ~ f^k.k', where 
2r7i£'k/(nsfc^) and 



This finally yields 



V TO V 2m ' 



Hh= ^kaj^flk- (10) 

k(fc<fco) 



Equation (9) has the same form as the Bogoliubov spectrum for a weakly coupled charged Bose gas, with the superfluid 
density rig taking the place of the condensate density. It will be used in the next Section as the dispersion relation of 
single-particle excitations in the evaluation of the phase-phase correlations from the single-particle Green's function. 
Equation (9) should also be contrasted with the dispersion relation of collective excitations, which within a weak- 
coupling theory is 



This relation tends to the 2D plasma frequency ujp = {2Trne^ /my^'^ at long wavelengths, in accord with the sum 
rule reported in Eq. (2). Notice that in the limit T 0, where we expect that Ug — > n, Eqs. (9) and (11) give the 
same result. According to the Gavoret-Nozieres theoremtj, single-particle and collective excitations coincide at long 
wavelengths in a Bose-condensed fluid at zero temperature. 

We conclude this Section by showing that the sum rule (3) is satisfied in our approach. The response function 
X^^t,^ (k, w) for the superfluid velocity is related to the phase-phase response by Xt,^t,^ (k, cj) — fc^X00(k, Li;)/m^. At 
long wavelengths we get from Eq. (8) 

lim xv^v^ (k, uj) = -^5^ ( 1 (12) 

k^o " " 2'm'^ilp ^ Qp + ie uj + \lp + is J 

where e = 0+ and flp = (27rnse^/TO)^/^. Insertion of Eq. (11) into the integral on the LHS of Eq. (3) leads immediately 
to the desired result. 
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III. ASYMPTOTIC BEHAVIOR OF THE SINGLE-PARTICLE GREEN'S FUNCTION 



($(k, U:M-k, ^ ——±rj;^ (14) 



We show in this Section that the one-body density matrix of the 2D-BCF at temperature T below a critical tem- 
perature Tc has a power-law decay, as a consequence of the correlations between phase fluctuations in the superfluid. 
We follow the method proposed for the neutral 2D gas in the work of PopovO (see also Fisher and Hohenber^) , 
which derives the power law by evaluating the single-particle Green's function at low momenta from the mean-square 
fluctuations of the phase $(r, r). 

More precisely, the phase fluctuations determine the single-particle Green's function G(r, r;ri,Ti) in the low- 
momentum regime (below the cut-off momentum fco) according to 

G(r, r; n, n) cx exp {-i ([$(r, r) - $(ri, ri)]^)}. (13) 

From Eq. (19.16) in Chapter 6 of Popov's bookH and using the dispersion relation for single-particle excitations given 
in Eq. (9) we find 

so that 

X I exp [i(k • r — cjt)] — exp [i(k • ri — LJTi)]p. (15) 

The summation over Matsubara frequencies in Eq. (15) can be carried out explicitly for ri = and for r = |r — ri| — > 
oo we find 

-([$(r,.)-$(r„n)n^^y coth (^14^ j y^^^ ^[1-Jo(M]. (16) 

Here, the temperature T is in units of e^/kB, L is a length scale of order 1/fco, and Jo{x) is the zero-order Bessel 
function. Finally, the expression for the mean-square phase fluctuation has the form 

\ ([$(r, r) - $(ri, ri)]2) ^ aln(r/L), (17) 



where the quantity a is given by 



^"^'''cothfil^iM^L (18) 



4 V n, / V Tr 



In conclusion, the one-body density matrix of the 2D-BCF decays to zero with the power law r~" as a consequence 
of the logarithmic correlations between phase fluctuations. In the limit of zero temperature the decay follows the law 
r~''=/4, which is the result found in the QMC study of Magro and Ceperleyi. 

Having reached their result from the finite-temperature formula given in Eq. (18), we can draw two main conclusions: 
(i) the off-diagonal order reported in the QMC work of Magro and Ceperleyo proves that the 2D-BCF is a superfluid 
at T = 0; and (ii) a quantum simulation study of the one-body density matrix at finite temperature would provide for 
the 2D-BCF an alternative approach to the superfluid density, to be compared with the methods that are currently 
available for its calculation (see e.g. the method based on the relation between the superfluid density and the mean 
square winding number0). Notice also that, as discussed in the book of Forsterlll, the equivalence between the 
definition of superfluid density that is provided by the single-particle excitation spectrum and the definition that may 
be extracted from transverse current correlations and rotation experiments remains a question of continued interest. 

As a final point, we should remark that our calculation has not included the possibility of vortex formation in 
the superfluid as the temperature is raised towards the critical temperature for superfluidity. Equation (18) should 
therefore be expected to become invalid near Tc. 
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IV. ALTERNATIVE APPROACH TO THE ONE-BODY DENSITY MATRIX 



We present in this Section an alternative calculation of the asymptotic behavior of the one-body density matrix, 
following an approach which is closely akin to the rcsummation method adopted in the book of PopovS and in the 
work by Kagan and coworkers! for neutral 2D Bose fluids in a quasicondensate state. 

We start from the definition p(r,r') — ^'(/;^(r)-0(r')) and use the expression of the field operators in terms of the 
hydrodynamic density and phase fluctuation operators to find 



p(r, r') =n,-\^b{v- r') + ^ (?7(r)'7(r')) + ($(r)$(r')) , 
to lowest order in the density and phase fluctuations. The transformations carried out in Section II yield 



($(r)$(r'))- ^ ak 

k(fc<fco) 



exp [ik • (r — r')] 



(19) 



(20) 



and a similar expression for {ri{r)r]{r')) , with f{E]s_) = (oIxl]^ ~ [exp (i^k/fes?") — 1] ^. The final result is 

1 



p{r) 



E 

k(fc<fco) 



471s ak 



n^ak 



/(i^k) 



exp [ik • r] 



(21) 



where r = |r — r'|. This expression has the same form as in the standard Bogoliubov approach, except that the 
superfluid density Us replaces the condensate density n^. Equation (13) can be re-obtained by an approximate 
resummation of the phase fluctuations to infinite order. 
We write the density matrix in the formi'i 



p(r) = p{r) exp [-A(r)] 

where A(r) collects all the terms that are responsible for the slow asymptotic decay of p{r). We find 

I/O . ^ 



(22) 



4 



1/2 



dx 



1 - Jq{xR) 



xg{x) 



1 



exp[2{ns/ny/^g{x)/{rsT)]-l 



(23) 



from Eq. (21), where R = 2[n,/ {nr'^J]^/'^{r /Iq) and g{x) = (1 + x^)^/^. We also find 



p{r) 



1 -r. 



(27r)2 



[1 — cos(k • r)] 



2niE^ 
xdx[l — Jf){xR)\ 



x^ 



(24) 



. 2g{x) g{x){exp [2{ns/ny/^gix)/{rsT)] - 1) 

In these equations the temperature is scaled in units of e^/kB- 

A power-law decay of p(r) can now be demonstrated analytically from Eqs. (22)- (24). The function Jq{xR) in the 
integrand in Eq. (23) provides a lower limit of integration going as r^^, while the upper limit is set by the cut-off 
wave number kg w 1/L for the applicability of the hydrodynamic Hamiltonian. The asymptotic calculation of the 
integrals yields 



p(r) ^ po(r/L)-" 

where the value of the exponent a is given by Eq. (18) and we have defined 

Pq = lim p{r). 



(25) 



(26) 



In the numerical calculations that are reported in the next Section we find that the ratio po/n becomes larger than 
unity at very low temperature, so that this quantity cannot be interpreted as a quasicondensate fraction. 
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V. NUMERICAL RESULTS 



In this Section we report some illustrative calculations of the one-body density matrix of the 2D-BCF, based on 
Eqs. (22) - (24). We need for this purpose to first evaluate the superfluid density as a function of temperature, for 
which we adopt the Landau theory based on damping of superfluid flow by emission of collective excitations. That 

isBi, 

n [exp(/3wk) - 1]^ 

with /3 = l/{kBT) and uj^ given by Eq. (11). This yields 

_ 2 x3exp[2g(x)/(r,r)] 

TJo {exp [2g(x)/(r,T)] - ip ^ ' 

with T in units of /ks- The results for the superfluid fraction ns{T)/n are shown in Figure 1 for some values of 
in the weak-coupling regime. 

It is worth noting that, as a consequence of the behavior shown by the superfluid fraction in Figure 1, the single- 
particle excitation energy i5k given in Eq. (9) starts for T = at the collective excitation energy Wk given in Eq. (11) 
and decreases with increasing temperature till its leading long-walength term vanishes at the critical temperature 
Tc- The plasmon gap in the collective excitation spectrum of Eq. (11) remains instead constant with increasing 
temperature. This behavior of the plasmon gap in the single-particle spectrum of a charged boson fluid has previously 
been found in the 3D case within the Hartree-Fock-Bogoliubov approximationEl, with small deviations from the 
predictions of the Gavoret-Nozieres theorem being present at T = 0. 

Our results for the one-body density matrix are reported in Figures 2 and 3, first at T = for several values of 
rs (Figure 2) and then at =0.1 for several values of T (Figure 3). In each of these Figures the right-hand panel 
shows a logarithmic plot of the numerical results for p{r)/n and a linear fit based on Eq. (25), using the value of 
the exponent a given by Eq. (18). The length L in Eq. (25), which is not determined by the asymptotic calculation, 
appears in a logarithmic scale as an additive constant and has been obtained from the linear fitting procedure at large 
r. Evidently the power-law prediction is in excellent agreement with the numerical results over an extended range of 
values for the reduced distance . 



VI. SUMMARY 



In summary, we have studied superfluidity in a weakly interacting 2D fluid of charged bosons with ln(r) interactions 
as a function of temperature, using a hydrodynamic reduction of the Hamiltonian that has allowed us to treat the 
system as a superfluid in the absence of a Bose-Einstein condensate at zero temperature. In the Appendix we show 
that the absence of a condensate at T = also follows for this system from the standard Bogoliubov approach. 

We have shown that the assumption of complete superfluidity for this system at T = accounts quantitatively for 
the QMC data of Magro and CeperleyS on the power-law decay of the one-body density matrix. In turn, computer 
studies at finite temperature would allow one to determine the behavior of the superfluid fraction and provide a full 
quantitative test of our theory. 
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APPENDIX. ABSENCE OF BOSE-EINSTEIN CONDENSATION IN THE 2D-BCF 



We show here that within the BogoHubov approximatioiJli there is no Bose-Einstein condensate in the 2D-BCF both 
at zero and at finite temperature. In contrast, in the 2D fluid of charged bosons with e^/r interactions a condensate 
is present at T = and the notion of quasicondcnsation can be used at T ^ olll. 

If we assume that the condensate fraction no/n is non-zero, then it is easy to show that the expression (9) for the 
single-particle excitation energy becomes 



2m 



(A.l) 



in the Bogoliubov approach. The corresponding equation for the condensate fraction reads 

f 1 _ v^r = ^ r^ixf 1 1 ] (A2) 

V nJKnJ 2 Jo 1 2x^g{x) ^ cxp [2{m/ny/^g{x)/{rsT)] - 1 J ^ ' ' 

where h{x) — 1 + 2x'^. 

The first term in the square brackets on the RHS of Eq. (A. 2) gives a contribution of order x~^ to the integrand 
as a; ^ 0, and hence the integral diverges. We conclude that no/n must vanish at any temperature, in agreement 
with the general sum-rule argument given by Magro and CeperleyEI. The property riQ/n = is consistent, on the 
other hand, with the Bogoliubov - des Gennes equations for the Bogoliubov amplitudes Wk and Wk- In this case these 
equations yield Wk = and Wk = 1, so that the particle density is related to the (now finite) chemical potential by 

-ll . (A.3) 



E 



exp 
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FIG. 1. The supcrfluid fraction n^/n as a function of reduced temperature T/Tc for four values of the couphng strength Vg- 
In absolute units the critical temperature takes the values 136.9, 41.3, 8.9 and 2.9 e^/fcs for going from 0.1 to 1. 



£ 0.7 





10 15 20 

r/lo r/lo 

FIG. 2. The one-body density matrix p{r)/n as a function of distance r (in units of the reference length Iq) in the 2D-BGF. 

Left: at T = for ra = 0.1, 0.2, 0.5, and 1 (from top to bottom). Right: on a logarithmic scale, the same numerical results 
are shown as dots while the predictions from the Green's function approach are shown as continuous lines. 
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FIG. 3. The one-body density matrix p(j)ln as a function of distance r (in units of lo) in the 2D-BCF. Left: at Ta = 0.1 for 
T — 0, 10, 20, and 40 e^/ks (from top to bottom). Right: on a logarithmic scale, the same numerical results are shown as 
dots while the predictions from the Green's function approach are shown as continuous lines. 



